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■ 1 INTRODUCTION 



ABSTRACT 

The aim of the present paper is to investigate a possible contribution of the rotation- 
powered pulsars and pulsar wind nebulae to the population of ultraluminous X-ray 
sources (ULXs). We first develop an analytical model for the evolution of the distri- 
bution function of pulsars over the spin period and find both the steady-state and the 
time-dependent solutions. Using the recent results on the X-ray efficiency dependence 
on pulsar characteristic age, we then compute the X-ray luminosity function (XLF) 
of rotation-powered pulsars. In a general case it has a broken power-law shape with 
a high luminosity cutoff, which depends on the distributions of the birth spin period 
and the magnetic field. 

Using the observed XLF of sources in the nearby galaxies and the condition that 
the pulsar XLF does not exceed that, we find the allowed region for the parameters 
describing the birth period distribution. We find that the mean pulsar period should be 
greater than 10-40 ms. These results are consistent with the constraints obtained from 
the X-ray luminosity of core-collapse supernovae. We estimate that the contribution 
of the rotation-powered pulsars to the ULX population is at a level exceeding 3 per 
cent. For a wide birth period distribution, this fraction grows with luminosity and 
above 10 40 erg s _1 pulsars can dominate the ULX population. 

Key words: methods: statistical - pulsars: general - stars: luminosity function, mass 
function - stars: neutron - X-ray: galaxies 



Ultraluminous X-ray sources (ULXs) are non-nuclear, point- 
like objects with apparent X-ray luminosity exceeding 
the Eddington lim it for a stellar mass black hole (see 
iFeng fc Sorial201ll . for a review). These objects were discov- 
ered by_the_Einsiem_s^tejhte in ne a rby star- forming galax- 
ies llLong fc van Speybroeckl Il983l; iFabbiand Il989l , Il988t 



iFabbiano fc Trinchieril Il987l ; IStocke et al.l Il99ll ). Observa- 
tions with Chandra and XMM-Newton satellites have ex- 
tended the sample of probable ULXs t o about 500 sources 
|Swartz et al.ll201ll ; IWalton et al.ll201ll ). 

There are several hypotheses about the nature of 
ULXs. The most popular models at this moment involve 
stellar-mass objects similar to SS 433 with the supercrit- 
ical regime of accretion a nd mild beaming with beaming 
factor 1/b = 4n/Q < -, 10 l|King et all l200ll ; iFabrikal 12004 
iPoutanen et alj|2007t ). or the accreting intermediate mass 



black holes (IMBH) with masses M 
IColbert fc Mushotzky|[l999l ). 
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^10 3 -10 5 M© (e.g. 



Many ULXs show spectral variability 

ijKaiava fc PoutanerJ I2009T ) typical for the accreting black 
holes. The presence o f soft thermal excesses sometimes seen 
in the ULX spectra (|Kaaret et all 120031 : [Miller et alj 12003! ') 
can be used as an argument of a large emission region size, 
which is either a signature of an IMBH or, alternatively, 
a large extended photosphere in a strong outflow from 
stellar-mass objects a ccreting at super-Eddington rates 
(|Poutanen et alj |2007| L The best IMBH candidates, the 
brightest ULXs, M82 X-l and ESO 243-49 HLX-1, show 
spectral states similar to those seen in Galactic sources 
llGladstone et al.ll2009l : IFeng fc Kaaretll2010l ; IServillat et all 
I2011IL but at higher luminosities. However, IMBHs cannot 
dominate the ULX population, because ma ny ULXs are 
assoc iated with the star-forming regions (|Swartz et al.l 
2009) and young stellar cl usters, but are clearly displaced 
from them by 100-300 pc (Zezas et al.ll2002l : iKaaret et all 
l2004lPtak et al.ll2006l : IPoutanen et all 12012] '). This in turn 
strongly argues in favour of the young, massive X-ray 
binaries as the ULX hosts that have been ejected from the 
stellar clusters by gravitational interactions during cluster 
formation and/or due to the supernova (SN) explosions. 
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It is very likely, however, that the ULX class is not 
homogeneous, but contains different kinds of objects. 

For example, some of the bright, steady ULXs could be 
young, luminous rotation-powered pulsars. Earlie r studies 
i|Seward fc Wand 1 19881 ; iBecker fc Truemperl Il997f ) suggest 
that X-ray luminosity of the pulsars is correlated with the 
rotation energy losses L = rjE I0 t. The efficiency rj, which 
defines the amount of rotational energy losses converted to 
the X-ray radiation, was found to be nearly constant. Later, 
using a more complete samp le of X-ray rotation-powered 
pulsars, IPossenti et alj (|2002l ) showed that the X-ray lumi- 
nosity depends on the rotational energy loss as a power law 



iPerna fc SteUal (2004) performed first Monte Carlo sim- 
ulations of the X-ray luminosity function (XLF) of rotation- 
powered pulsars. In order to describe the luminosity evo- 
lution of the pulsars together with the evolution of the 
spin period due to the magnetic-dipole radiation losses, 
they used the effic ie ncy - characteristic age dependence 
from IPossenti et"al] (|20()2h . They considered the distribu- 
tion functions of pulsars o ver the magnetic fi e ld an d the 
birth spin period given by lArzoumanian et al.l l|2002T l and 
showed that rotation-powered pulsars can be very bright X- 
ray sources with luminosities L > 10 39 erg s _1 . 

Recent investigation of the X-ray pr operties of 
rotati onal-powered pulsars conducted by IVink et al.l 
(2011) revealed a more complicated efficiency-age de- 
pendence. Using the new d ata from Chandra observatory 
l|Kargaltsev fc Pavlovl [2008h . they find that radiative effi- 
ciency is not constant for pulsars with age < 1.7 x 10 4 yr, 
but depends on the characteristic age. These new results 
may strongly affect the XLF, increasing the number of the 
most luminous pulsars. 

In the present paper we develop a model for the XLF 
of the rotation-powered pulsars, taking into account the re- 
cently discovered efficiency-age dependence. In Section[3]we 
present the analytical model describing the evolution of the 
pulsar periods and find both the steady-state and the time- 
dependent solutions. Section is devoted to the observa- 
tional constraints on the model parameters for the birth 
period and magnetic field distribution that can be obtained 
from the core-collapsed SNe and the observed XLF of the 
sources in the nearby galaxies. In Section [4] we obtain the 
birth period and magnetic field distributions for the bright- 
est pulsars and estimate the possible contribution of young 
pulsars to the ULX population. We summarize in Section [5] 



2 MODEL 

2.1 Basic equations 

A pulsar is described by two parameters: its birth period 
po and the magnetic field B, which is assumed to be con- 
stant over its lifetime. We consider lognormal distributions 
for both B and po, with the probability density for the dec- 
imal logarithm log x in the following form: 



G{\ogx\\og{x),a x ) = 



log 2 W(*)) 
2<j2 



2ir a x 



(1) 



The mean and the standard deviation (scale) for the two 
distributions are (log(S), o\b) and (log(po), 0j>). 



The pulsar period at a given age is calculated using a 
simple model, where the rotational energy losses ar e dom- 
inate d by the magnetic dipole radiation (see e.g. iGhoshl 
l2007t) : 



En 



-Jfifi 



2-R p 2 4 



(2) 



where Q — 2-k/p is the pulsar rotational frequency, R is the 
neutron star radius and / is its moment of inertia. We ignore 
the factor depending on the angle between the dipole and 
the rotational axis to be consistent with previous studies. 
The evolution of the pulsar period and the frequency are 
described by equations 



■ 1 R 2 
pp = -olB , 

where 

16tt 2 R 6 



il = --^B 2 tf. 

Sit 2 



3 7c 3 ' 

The time-dependence of the pulsar period is then 



p(t) = sJpl + aBH, 



with the characteristic spindown age 

= 5 x 10 14 p 2 Br 2 2 s, 



c ~ 2p aB2 



(3) 



(4) 



(5) 



(6) 



where we assumed Ins = 1 and Re = 1 (i.e. a « 2 x 10 -39 
cgs) and used standard notations Q = 10 X Q X in cgs units. 
On the B-p plane, using equation © we can identify the 

2/2 

lines of constant characteristic age B12 = 3.9 ps r c , y r (see 
dotted lines in Fig. [!}. 

We estimate the X-ray luminosity L of a pulsar (and 
a pulsar wind nebu la, PWN) f r om it s period and period 
derivative following IVink et alj (|201ll ). For simplicity, we 
approximate the efficiency-age dependence with a simple 
relation: 



(7) 



{f?o, if r c < n, 

TIG (Tl/rc) 2 , if Tl < Tc < T2, 
10" 4 , ifr c >r 2 , 

where ti = 170 yr x ri^ 1 ^ 2 and T2 = 1.7 x 10 4 yr. This is 
equivalent to 

f 4 x 10 31 p" 4 Bi 2 2??o erg s" 1 , if r c < n, 



4 x 10 21 p" 8 Bf 2 erg s" 1 , if n < t c < r 2 , 



(8) 



4 x 10 2 'p- #12 erg 



if T c ^ T2 . 



The value of the maximal efficiency r/o is not well de- 
fined because of the lack of young pulsars in the Milky 
Way. The data seem to indicate that it is at least 0.3 
l|Kargaltsev fc Pavlovl20"0sl : IVink et al . 2011), which we take 
as a lower limit. In principle, it can even exceed unity, be- 
cause of the beaming of the pulsar radiation. 

In the B-p plane, we can identify the lines of constant 
luminosity (see solid lines in Fig.[l}. Depending on the range 
of r c , these lines have different slopes (see equation JHJ): 



Bio = 



5 x 10~ 3 V 1/2 , if T C < Ti or T C > T 2 , 



0.08 p 



4/3 r l/6 



if Ti < T c < T 2 . 



(9) 



The line L — const intersects with the line r c = 
ti at a point 1 with coordinates (p-3,1, _Bi2,i) = 
(60L^ 9 1//2 77q' /4 , l8L~g^ 2 rio), while an intersection with the 
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Figure 1. Contour plots of constant X-ray luminosity log L (solid 
lines) and the characteristic age r c (dotted lines) on the plane 
magnetic field - period for tjq = 1. 



line 



T2 occurs at point 2 (p_ 



,B 



12,2 ) 



(60L 3 - 1/2 ,1.8L-3 1/2 ). 

As the pulsar period increases, its luminosity drops. If 
B12 > -Bi2,i or B\2 < B\2,2, a pulsar crosses the line of 
a given luminosity being at the constant efficiency branch 
r\ = ?)o or r\ — 10 -4 , while for Bi2,2 < B12 < -812,1 it occurs 
at the decaying branch of 77. Thus for the fixed magnetic 
field, the pulsar period at a given luminosity L is 



1A ni/2 r -i/* !/ 4 
14 B X 2 L 39 Vo , 

,. = { 6.7 BH 4 L; g 1/s , 

45 B\fL-V\ 



if -B12 > -812,1 , 

if Bl2,2 < B12 < Bl2,l, 
if B±2 < £>12,2 ■ 



(10) 



2.2 Steady-state distributions and the differential 
luminosity function 

2.2.1 Steady-state period distribution 

The evolution of the distribution function of pulsars over the 
period N(p) = dN/dp (for a given magnetic field B) can be 
described by the following evolution equation: 



(11) 



with the source function describing the production of new 
pulsars per unit period and time given by 



Q(p) = N 



1 



pin 10 



G(logp;log(po),crp) 



(12) 



and N — j Q(p)dp is the total production rate per unit 
time. Equation (|11[) can be solved analytically and in the 
steady-state the solution takes the form: 
1 Cf 



N( P ) = - I Q(p'W- 
V Jo 



(13) 
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Figure 2. Differential XLF of rotation-powered pulsars for the 
pulsar birth rate N yI = 0.01. Dashed red line shows the XLF for 
the fixed magnetic field B12 = 4 and very small birth periods with 
log(po) = —2.5 and a v = 0.2. The maximum efficiency is assumed 
to be 7/0 = 1. The dash-dotted green line corresponds to the larger 
birth periods with log(po) = —2, and the XLF for even larger 
log(po) = —1.5 is shown by the dotted line. The XLFs averaged 
over the magnetic field distribution with parameters log(B) = 
12.6 and ag = 0.4 for the birth period distribution parameters 
log(po) = —1.7 and a p = 0.2 are presented by the black solid 
and dashed lines for the maximum efficiency of rjo = 1 and 0.3, 
respectively. 



It reduces to 



p 

aB 2 



1 + erf 



log(p/(P(>)) 
<T p v2 



(14) 



for the source function given by equation (|12[) . At periods 
which are much larger than the initial periods we get 



N(p) = — = 3.2 x 10 7 Br 2 2 N yl p, 
P 



(15) 



where N yT is the pulsar birthrate per year. 

If the magnetic field and the birth period distributions 
of the pulsars are lognormal, the steady-state period distri- 
bution averaged over the magnetic field distribution is given 
by 



{N(p)) B = N 



a(B) 1 



-P 



1 + crf( lo g (p/(po)A 
a p V2 J 



(16) 



For the periods p much larger than (po), the distribution has 
a form: 

.2 .2, 



(N(p))i 



3.2 x 10 



7 e 



„2ctr In' 1 10 



<Bi2> 2 



■ Wyr P, 



(17) 



where we used the relation (B^) = (B) e exp[(Ca s lnlO) 2 /2]. 



2.2.2 Luminosity distribution 

For a given magnetic field, the XLF can be obtained from 
the period distribution function via transformation 



L N(L) = pN(p) 



d log p 



dlog L 



(18) 
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For the power-law dependence of luminosity on period L — 
C p~ 7 , we get 

LN(L) = 3.2 x 10 7 Br 2 2 Nyr -{L/C)~ 2h . (19) 

7 

For the constant X-ray efficiency rj (i.e. 7 = 4), we then 
easily get from equation 



LN(L) = 5 x 10 3 B^V 1/2 N yr L~g /2 . 



(20) 



Thus for very young, rapidly rotating, luminous pulsars as 
well as for the old pulsars the distribution will follow that 
law (see Fig. [5} In the intermediate regime for n < r c < Ti, 
7 = 8 and the XLF follows a shallower dependence: 



LN(L) = 3.2 x 10 2 B^ /2 N yt L 3g 



1/4 



(21) 



According to equation @, the high- luminosity break is ex- 
pected at L — 4 x 10 41 B^ 2 erg s _1 . Smaller mean periods 
and larger magnetic fields lead to a larger initial luminosity 
and therefore to a larger number of luminous sources. A cor- 
responding low-luminosity break is at 4 x 10 33 -B-^ 2 erg s _1 . 

The examples of the XLF normalized to the pulsar birth 
rate N yT are presented in Fig. [2] The XLF has a complex 
shape reflecting the behaviour of the X-ray radiative effi- 
ciency. The XLF for the fixed B has sharp features reflect- 
ing breaks in the derivative of the function L(p) given by 
equation ©. These breaks are unlikely to be observed, be- 
cause the actual efficiency-age dependence (|7J is likely to be 
smooth. 

The luminosity function of pulsars with the magnetic 
field distribution G(log B; log(B), erg) can be obtained by 
averaging the XLFs over that distribution. In this case, the 
sharp features also disappear (see solid line in Fig. [2J). In 
the range of luminosities corresponding to the constant effi- 
ciency, we get 



L N(L) = 5 x 10 3 



s (<r B lnlO) 2 /2 

(B12) 



1/2,7 r-1/2 
77 JV yr _l 38 



(22) 



In the intermediate range of luminosities, the power-law seg- 
ment has the following form: 



LN(L) = 3.2 x 10 



(agInlO) 2 /8 

Bi2 y /2 "yr ^38 ■ 



(23) 



If the birth period distribution has a peak at rather large 
periods, the XLF has a cutoff before the high-luminosity 
power-law segment actually starts (e.g. see black solid line 
in Fig. [2]| . Decreasing the maximum efficiency 770 leads to a 
smaller cutoff luminosity (compare solid and dashed black 
curves in Fig. while the intermediate power-law barely 
changes. 

Radiation from a pulsar may be confined w ithin a 
narrow beam ~ 1 str ijTauris fc Manchester! 1 19981 ) corre- 
sponding to the beaming factor b — £1/4tt = 0.1. How- 
ever, young pulsars have PWN, which are more isotropic. 
The ratio of the observed luminosities of the nebula to 
the pulsar has a large spread, but t ypically is of the or- 
der unity l|Kargaltsev fc Pavlovll2008h . This argues against 
strong beaming and therefore we take b > 0.3. The normal- 
ization of the observed luminosity function scales linearly 
with the beaming factor. 



2.3 Non-stationary solution 

2.3.1 Evolution of the period distribution 

In order to determine the distribution of pulsars over the 
period at a given pulsar age, we need to solve the time- 
dependent evolution equation (JTTJ : 



dN( P ,t) + d 



in 



9- p [pN(p,t) 







(24) 



with the following initial condition at zero-age: 

N(p,t = 0) = N o (p). (25) 

According to equation ((5JI, the period derivative can be ex- 
pressed as p = aB 2 /2p. Equation (|24fl conserves the total 
number of pulsars N. Its solution is 



N(p,t)=p 



N {po) 
Po 



(26) 



where po = \J p' 1 — atB 2 . The solution is only defined for 
p ^ BVed. (27) 
Function N(p, t) has the mean 



(p)(t) = J P N{p,t)dp = j yfpJ+rtB* N Q (po)dp , (28) 


and the variance 

a 2 p (t) = J (p-(p)) 2 N(p,t) dp 


00 „ 
= J (jpl + atB 2 - (pfj iVo(po) d Po . (29) 



At large ages, the mean becomes 

<p3> 



<p)(t)My/atB + 

2V can 

and the variance is 

Thus, the solution becomes the delta-function: 



N(p, t — > 00) = TV 5{p — VatB) 



(30) 



(31) 



(32) 



This can be easily understood from equation ([5|, which de- 
scribes the evolution of the pulsar period with time. At large 
age, the time-dependent term becomes much greater than 
the value of the initial spin period. Thus, every pulsar with 
a given magnetic field at a given age has the same period 
pit) = \fatB. Therefore, at large ages the period distribu- 
tion of the pulsars does not contain any information about 
the initial one. The characteristic timescale at which the 
information about the initial distribution is lost can be es- 
timated as: 



t ' 



poV 
B) 



1600 



Po, 



B1 



yr. 



(33) 



where po,-2 is the birth period expressed in 10 ms. 

We can now obtain the solution averaged over the mag- 
netic field distribution 
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Figure 3. Evolution of the (normalized) period distribution with 
time. The black solid lines show the evolution of period distribu- 
tion with initial o p = 0.2 averaged over magnetic field distribution 
with o B = 0.4 for r = 0, 1, 10, 100 (from left to right). The blue 
dotted line shows the period distribution for a specific magnetic 
field (/3 = 1) at r=l, 10, and 100. The distribution becomes very 
narrow at large r, peaking at p = y/r. The red dashed lines rep- 
resent the asymptotic solution H37H at large r, which just reflects 
the lognormal distribution of the magnetic field. 



oo 

{N(p,t)) B = J N(p,t) G(logB;log(£), a B ) dlogB. 



It is useful to introduce dimensionless variables 



(34) 



t. 



(35) 



V a B ((BY 

and find the solution as a function of dimensionless period 
p and time r, such as (N(p,t))b = {N(p,t))s X (po) and 
/ {N(p, t)) B<lp = N. For the lognormal distribution of both 
magnetic field and birth periods, we get 

P 



{N( P ,t))b=N 



2ira p a B In 10 



(36) 



dp 



P(P 2 ~ rP 2 ) 



exp 



log y p 2 - w log 2 /? 



2a2 



7a% 



An asymptote at large r can be easy obtained directly sub- 
stituting (|32[1 to equation (I34[) : 

N 



{N(p,T^oo))i 



-G(logp;log y/T, a B ) 



(37) 



pin 10 

This implies that the mean of the distribution pN(p) in- 
creases with time as (p) — y/r and the dispersion is com- 
pletely determined by the width of the magnetic field distri- 
bution gb- Evolution of the period distribution is presented 
in Fig. [3] If a p < aB, then at large r the period distribu- 
tion is wider than the initial one, while in the opposite case, 
a v > as, the period distribution becomes narrower. 

2.3.2 Evolution of the luminosity distribution 

In order to obtain the luminosity distribution, we take the 
time-dependent solution for the period distribution (I26[) . use 



Figure 4. Normalized luminosity distribution at birth for differ- 
ent values of birth period with <r p = 0.2, log(_B) = 12.6, erg = 0.4 
and r?o = 1 (black lines). Blue dotted and red dashed lines show 
the luminosity distribution with gb = 0.2 and a p = 0.4, respec- 
tively, for log(po) = — 1.5 (other parameters as above). 
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Figure 5. Evolution of the pulsar normalized luminosity distri- 
bution for the initial distribution with the following parameters 



log(po) 



-1.5, o-p = 0.2, log(B) = 12.6, o B = 0.4 and % = 1 



(black solid lines). Blue dotted and red dashed lines shows the 
luminosity distribution at 1000 yr for erg = 0.2 and o p = 0.4, 
respectively (other parameters as above). The luminosity distri- 
bution becomes more symmetric at large age. 



the transformation (|18[) and average the derived expression 
over the magnetic field. Resulting distribution and its evolu- 
tion is presented in Figs[4]and[5] respectively. As it is clearly 
seen from Fig. [4] the initial luminosity distribution of the 
pulsars can be multimodal. Every mode of the distribution 
is related to the different regime of the luminosity-period de- 
pendence. Also, the luminosity distribution at birth may re- 
veal the narrow spikes, related to the breaks in the dp(L) /dL 
derivative. The luminosity distribution is broader for larger 
cr p and aB- The distribution becomes narrower and more 
symmetric as the time increases (Fig. [5]l . This happens be- 
cause at birth, pulsars can operate in the different regimes 
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Table 1. Birthrates and parameters of the magnetic field and the 
birth period distributions for rotation-powered pulsars. 



# 


log(po) 


Op 


log(B) 




N a 


Reference 6 


1 


-2.3 


0.3 C 


12.35 


0.4 


0.0013 


1 


2 


-1.52 


0.0 


12.75 d 


0.33 d 


0.01 e 


2 


3 


-0.52/ 


0.8-f 


12.65 


0.55 


0.028» 


3 


4 


-1.7' 1 


O.l' 1 


12.6 


0.1 


0.01 


4 



a Birth rate of p ulsars in the Milky Way pe r year 



b Refe renc es: (1) [A rzoumaniaji_et_al 



200211 : (31 IFaucher-Giguere , 
201ll l. 



Kaspi (2006); (4) 



(2002); (2) 



Gonthi er et al.l 
TakateTeTlrt] 



c Ref. 1 gives the lower limit on cr p of 0.2. Taking a broader 
distribution with cr p > 0.3 does not affect the results. 
d Parameters for the lognormal distribution were estimated by 
fitting a more complex distribution adopted in ref. 2, see their 
Table 1 and eq. (1). 

e Value from the first line of Table 8 of ref. 2. 

f Parameters for the lognormal distribution were estimated by 

fitting a Gaussian distribution adopted in ref. 3, see their Table 

8. 

9 Birthrate from Table 8 of ref. 3. 

" The lognormal distribution approximates the flat distribution 
in the 20—30 ms range adopted in ref. 4. 



with a Gaussian with the mean po~0.1 s and width <r~0.1 s. 
However, this result is inconclusive, because the number of 
objects in the used sample is not large enough to derive the 
exact shape of the period distribution. 

The analysis of the observed lum inosity distr i bution 
of the historical core-collapse SNe by iPerna et alj (|2008l ) 
showed that th e pred icted number of bright pulsars in the 
IPerna fc Steilal l|2004 ) model is much larger than the ob- 
served number of luminous SNe. This discrepancy is re- 
lated tc^tlie_assjuTi«d_very __short (5 ms) mean birth period 
from lArzoumanian et all l|2002t) . On the other hand, using 
parameters of the pulsar magn e tic field distributio n from 
IFaucher-Giguere fc Kaspl (|200fj| ), IPerna et all (|2008h found 
that in order to satisfy the observed luminosity distribution 
of the historical core-collapse SNe, the birth period of the 
pulsars should be larger than 40-50 ms. 

In the f ollowi ng sections we repeat the analysis by 
IPerna et al. (2008) using a different efficiency-age depen- 
dence given by equation (J7J) as well as using the new 
data that became available after 2008. We also obtain con- 
straints on the pulsar birth period distribution by compar- 
ing the simulated pulsar XLF with the observe d XLF of the 
bright sources in the nearby galaxies derived bv lMineo et al.l 
(|2012T). 



of conversion of the rotational energy losses to the X-ray 
radiation, depending on the spin period and the magnetic 
field distributions. With time, all pulsars move towards the 
same regime, where the efficiency is constant ~ 10 -4 (see 
equation (0). 



3 OBSERVATIONAL CONSTRAINTS ON 
MODEL PARAMETERS 

3.1 Previous determination of magnetic field and 
birth period distributions 

Distributions of the pulsars over the magnetic field and the 
birth period were investigated in several papers based on 
the analysis of the observed rad io (|Arzoumanian et "afll2002i ; 
IFaucher-Giguere fc Kaspll2006l) an d the gamma-ray pulsars 
( Gonthier et al. , 2002; Ta kata et al.ll201ll ). Parameters of the 
magnetic field distribution are similar in all these studies 
lying in the range log(B) = 12. 35-12. 75, <t s =0. 1-0.55 (see 
Table[l}. However, parameters of the birth period distribu- 
tions are significantly different: the mean logarithm log(po) 
varies from —2.3 to —0.5 (i.e. periods in the range from 5 
to 200 ms) and the width a p varies in the range 0-0.8 (see 
Table[T] the parameters were estimated by fitting the lognor- 
mal distribution to the actual distributions adopted by the 
authors). This difference in the birth period distributions 
is most likely caused by a low sensitivity of the considered 
models to the birth period. As it was shown on Section [21 
the period distribution of the pulsars at large time does not 
contain information about the birth periods. Therefore, in 
order to determine these parameters we have to use only 

young pulsars. 

Recently, IPqpov fc Turollal (|2012T l have presented new 
estimates of the birth periods based on a sample of radio 
pulsars associated with the SN remnants. They showed that 
the distribution has to be rather wide, and it is consistent 



3.2 Constraints from the luminosity distribution 
of core-collapse SNe 

IPerna et al.l l|2008l ) proposed that constraints on the birth 
period distribution can be obtained by comparing the ob- 
served luminosity distribution of historical core-collapse SNe 
with the simulated pulsar XLF (for the given ages), consid- 
ering that the most probable remnant of the core-collapse 
SN explosion is a neutron star. 

One of the important questions is the earliest age at 
which SNe can be used to derive the observational luminos- 
ity distribution that would reflect the XLF of the bright- 
est rotation-powered pulsars. There are two main issues 
here. The first problem is the high optical depth of the 
S NR shell at the earliest stages of its expansion. According 
to IChevalier fc Franssonl (|l994l ). the optical depth of SNR 
changes with time as 

(38) 

for typical parameters of the SN explosion (kinetic energy 
Eo ~ 10 erg, mass of the ejecta M e j ~ 10 Mq , and typical 
photon energy E ~ 10 keV). Therefore, the S NR shell be- 
comes optically thin in about 3 years. However. IPerna et alj 
(2008) limited their SNe subsamples by the minimal age of 
10 and 30 yr. As a result, they did not include the most 
luminous SNe in their analysis. Furthermore, because of the 
finite size and rapid expansion of the shell, the diffusion time 
of photons in the SNR shell can be small enough to make 
the X-ray radiation of the central pulsar visible at even ear- 
lier times. Considering that the number of scatterings in the 
SNR shell of optical depth r ^> 1 scales as N ~ r 2 , we can 
estimate the escape time of the photons from the shell: 

t (t)~ T{t)RsMt \ (39) 
c 

where J?sNii(t) is the size of SNR. Because of the free ex- 
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Table 2. X-ray luminosities of historical SNe 



Table 3. Upper limits fo r the X-ray luminosities of historical SNe 
(from iPerna et af1l2008h . 



SN 



Age (yr) a logL 



References" 



1979C 

1986E 

1986J 

1988Z 

1990U 

19941 

1995N 

1996cr 

1998S 

1998bw 

1999ec 

2001em 

2001gd 

2001ig 

2004C 

2005ip 

2005kd 

2006jd 

2008ij 



26.8 
19.6 
21.2 
15.5 
10.9 
8.2 
8.9 
4.2 
3.6 
3.5 
5.9 
4.7 
1.1 
0.5 
3.1 
1.3 
1.2 
1.1 
0.56 



on ,, + 0.06 

oo,,o -0.07 

38.15+.^ 

QO QQ~l _ 0-02 

oo.yo_ Q 03 
39.46lg^ 

39.63t°;°? 
39.28+. 

38.60+|i° 

39 4O+0-05 

40 7fi+ ' 08 

39.oolg;lJ 

37 54+0-20 

38.oot°;JJ 

40 20+ ' 07 
41-411 

39.001°;^ 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
1 
3,1 
5 



a Ages of SNe were calculated from the detection times listed at 
the website of the IAU Central Bureau for Astrono mical 

Telegrams, excep t for SNe fromlPerna et a h| d2008|). 

b Referenc es: (lllPerna et al.l [(2008T ) : (2lllmmler fc Pooled 



2007 1; (3) 
mB); (5) 



Immlcr et al. (2007); (4) Dwarkadas ■ 
Immler et ail (2003). 



Gruszki 



lev 
zko 



pansion stage for SNe with ages < 
will increase with time as 



100 yr, the size of SNR 



(40) 



Therefore, the diffusion time is 
O05 

U 



yr. 



(41) 



The diffusion time is equal to the age of SN at t ~ 0.2 yr 
and later it is always smaller. Therefore, a significant frac- 
tion of the SNR radiation may be produced by the central 
pulsar, because the luminosity of a typical newborn pulsar 
may achieve 10 40 -10 41 erg s — 1 , which is comparable to the 
highest observed SNe luminosities in our sample. 

The second issue is the fallback accretion onto a neu- 
tron star during early phases of SN explosion. According to 
IChevalierl (1 19891 ). radiation from the central pulsar begins to 
diffuse through the accreting matter when a reverse shock 
radius reaches the radiation trapping radius. It happens at 
> 0.5 yr after the SN. Therefore, we can expect that the 
central pulsars will contribute to the total luminosity after 
0.5-1.0 yr. This estimate is close to the limit coming from 
the diffusion time arguments. Thus, we will use the minimal 
age t = 0.5 yr. 

Another important question is a selection effect, which 
may strongly affect the observed luminosity distribution 
of SNe, because the younger is the source the brighter it 
is and the higher is the probability for it to be detected. 
For example, mo st of t he SNe luminosity measurements 
from IPerna et al l (|2008l ) are upper limits, because those 
sources are quite faint. Only 19 brightest and youngest 



SN 




log L 


1923A 


77.3 


35.78 


1926A 


75.3 


37.15 


1937A 


67.3 


37.11 


1937F 


62.1 


36.43 


1940A 


63.0 


37.00 


1940B 


62.6 


36.93 


1941A 


60.2 


36.74 


1948B 


55.1 


35.67 


1954A 


48.9 


35.20 


1959D 


41.6 


37.34 


1961V 


38.3 


37.79 


1962L 


41.2 


37.67 


1962M 


40.3 


35.57 


1965H 


37.7 


38.18 


1965L 


37.8 


36.76 


1968L 


32.0 


36.18 


1969B 


32.6 


36.58 


1969L 


30.3 


37.68 


1970G 


33.9 


36.69 


1972Q 


30.5 


38.48 


1972R 


31.9 


35.86 


1973R 


25.9 


37.89 


1976B 


26.2 


37.95 


1980K 


24.0 


36.81 


1982F 


22.6 


36.04 


1983E 


19.0 


37.66 


19831 


17.8 


36.23 


1983N 


16.8 


36.74 


1983V 


19.1 


37.85 


1985L 


14.9 


37.91 


19861 


17.1 


38.48 


1986L 


18.9 


38.15 


1987B 


14.1 


38.18 


1988A 


12.3 


37.38 


1991N 


11.8 


37.62 


1993J 


8.1 


38.00 


1994ak 


7.1 


37.57 


1996ae 


5.9 


37.79 


1996bu 


6.6 


37.32 


1997X 


6.1 


37.34 


1997bs 


2.5 


38.46 


1998T 


5.2 


38.30 


1999dn 


1,1 


37.77 


1999el 


5.6 


38.75 


1999em 


1.0 


37.15 


2000P 


7.2 


39.08 


2000bg 


1.3 


39.15 


2001ci 


2.5 


37.70 


2001du 


1.3 


37.58 


2002ap 


0.9 


36.49 


2002fjn 


4.7 


39.11 


2002M 


3.1 


38.88 


2003dh 


0.7 


40.70 


2005N 


0.5 


40.00 


2005at 


1.7 


38.48 


2005bf 


0.6 


39.78 


2005gl 


1.6 


39.53 
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Figure 6. (a) Cumulative normalized luminosity distributions for SNe (cither measurements or upper limits; pink histogram) with ages 
t > 0.5 yr. Black solid, red dash-dotted, green dashed and blue dotted lines correspond to average distributions for models 1—4 from 
Table [JJ respectively. Thick lines are for the case r)o = 1 and the thin lines are for rjo = 0.3. Here no beaming is assumed (b = 1). (b) 
Allowed region for the parameters of the birth period distribution. Parameters along the red and blue lines satisfy condition (1421 in 90 
and 68 per cent cases, respectively. The solid, dashed, dotted and dot-dashed lines correspond to different pairs of (i>, rjo) = (1,1), (1,0.3), 
(0.3,1) and (0.3,0.3), respectively. Regions to the right of these lines satisfy the data with higher probability. Calculations are performed 
for the average values of the magnetic field distribution log(_B) = 12.6 and erg = 0.4. Positions of parameters listed in Table[T]are marked 
by different symbols. 



sources have actual measurements of luminosity. Compi- 
lation of Dwarkadas fe Gruszkol (|2012h contains additional 
eleven sources with known luminosity, but there are only 
four sources with ages > 0.5 yr. In addition, because some 
fraction of the SNe X-ray luminosity is not related to the 
pulsar or PWN, also the actual X-ray detections here should 
be treated as upper limits on the pulsar luminosity. For the 
analysis we use the data on the ages and the X-ray lu- 
minosities o f core- collapse SNe (with ages > 0.5 yr) from 
IPerna et alj (|2008l ) with the addition of the new measure - 
ments from the compilation of Dwarka das fc Gruszkol (|2012| ) 
(see Tables [2] and [3} . The cumulative histogram of upper 
limits is shown in Fig. [6] by the bold pink line. 

We follow the recipe described in Section [2.3.21 and cal- 
culate the luminosity distribution of 76 pulsars for the ages 
of SNe listed in Tables [2] and [3] We then construct the 
average normalized cumulative distribution of luminosities 
/(< L) and compare it to the observed distribution. In the 
absence of beaming (i.e. b = 1), from Fig .[6]we see that only 
model 3 (|Faucher-Giguere fe K aspi 2006) satisfies the upper 
limit distribution, f or the maximum eff iciency 770 between 
0.3 and 1. Model 2 (|Gonthier et al.ll2002T l is also reasonably 
close, especially for rjo = 0.3. 

An additional effect appears if beaming is significant. 
Then most of the pulsars which appear to be faint in the 
X-rays, in reality could be very bright sources, but beamed 
away from us. Furthermore, about 10 per cent of SNe pro- 
duce a black hole after explosion instead of a neutron star 
l|Heger et al.ll2003l ), which can be modelled as an additional 
multiplicative beaming factor 0.9. The cumulative normal- 
ized luminosity function in that case would start from the 
value 1 — 0.96 at the low-luminosity end. For example, for a 
smaller beaming factor b = 0.3, most of the bright pulsars 
would be undetected. In that case, models 2 and 3 satisfy the 



upper limit distributions, mo del 4 is only marginally co nsis- 
tent with them, but model 1 IjArzoumanian et al. 2002) still 
contradicts the data. 

We can also find more general constraints on the pa- 
rameter set (log(po), (jp). With some high probability, the 
cumulative model distribution should be above the observed 
histogram of upper limits /(< L) at any luminosity. We can 
formalize this condition by computing the fraction of model 
distributions that satisfy the constraint 



mm 

L 



/pul S ar S (< IQ((logpo), (J v ) 
/SNc(< L) 



> 1. 



(42) 



Using Monte-Carlo method, we simulate 3000 sets of 76 pul- 
sars at given ages (see Tables [2] and [3} which follow given 
magnetic field and the initial period distributions, account- 
ing for beaming and for the 10 per cent black hole fraction. 
We then find the dependence cr p ((po)), which satisfies condi- 
tion (142[) in 90 and 68 per cent cases. The results depend on 
b and 770 (see Fig.[6p). The mean birth period of the pulsars 
cannot be shorter than ~ 15-30 ms for any a p . For larger 
dispersion a p ~ 1, the limiting value is between 60 ms and 
1 s, depending on the parameters and the sought probability. 
The constraints depend much stronger on beaming than the 
assumed maximal efficiency 770, because they come mostly 
from the low-luminosity SNe, whose number is not affected 
much by variation of rjo- Our constraint s on the mean per i- 
ods are consistent with those derived bv lPerna et all ( 20081. 
who found (po) >40-50 ms. 
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Figure 7. (a) Cumulative luminosity distribution of the pulsars normalized to the unit SFR and unit beaming 6 = 1. Black solid, red 
dash-dotted, green dashed and blue dotted lines correspond to models 1—4, respectively. The thick and thin lines correspond to the 
maximum efficiency of t)q = 1 and 0.3 , respectively. The pink solid straight line represents the average cumulative XLF of sources in 
the nearby galaxies l|Mineo et al.ll2012f ) taken in the form i144[) . (b) The predicted fraction (above a given luminosity) of pulsars in the 
high-luminosity end of the average XLF. Each line corresponds to the same models as in panel (a). 



3.3 Constraints from the XLF for sources in 
nearby galaxies 

3.3.1 Averaged XLF of nearby galaxies 

The averaged XLF of the bright sources in nearby star form - 
ing galaxies was recently obtained by [Mineo et al.l (|20 12f ) . 
The star formation rates (SFRs) in the galaxies of their 
sample are spread in a broad interval between ~0.1 and 
~ 100 Mq yr" 1 . The observed XLF is well approximated by 
a power- law: 

^- = 1.88 x Lgg 1 - 59 x SFR[M Q yr" 1 ]. (43) 

iMineo et all |2012t ) introduce the cutoff at L C ut,38 = 10 , be- 
cause of lack of statistics at luminosities above 10 41 erg s -1 . 
Here we do not introduce the cutoff and integrate the rela- 
tion (|43[> to infinity to derive the cumulative distribution: 

jV(> L 38 ) = 3.2 x Lag 0,59 x SFR. (44) 

3.3.2 Comparisons of the pulsar and observed XLF 

In order to make the comparisons between the pulsar 
XLF and the observed XLF of sources in the nearby 
galaxies, we first need to find the relation between the 
pulsar birthrate and the SFR. We assume the Galactic 
S FRmw = 2 Mo yr -1 , in ac cordance with the recent study 
of IChomiuk fc Povichl l|201ll ). However, using luminous ra- 
dio SN remnants and the X-ray point sources, these authors 
found that the Milky Way deviates from SFR expectations 
at the 1— 3cr level, hinting that the Galactic SFR is overesti- 
mated or extragalactic SFRs need to be revised upward. 

The estimations for the birth rate of pulsars in the 
Milky Way differs by an order of magnitude in various pa- 
pers (see Table [TJ), and on average is about 0.02 yr" 1 . The 
conversion between pulsar birthrate and the SFR can be 
expressed as follows : 



N [yr" 1 ] = k x SFR [M yr" 1 ] , (45) 

with the conversion factor k varying between 0.0007 (in 
model 1) and 0.014 (in model 3), with the mean of about 
0.01. 

Using the conversion formula (|45[1 we can now produce 
the cumulative luminosity distribution of pulsars normal- 
ized by the SFR a nd compare it to the observed XLF from 
IMineo et all (|2012T l. The XLFs calculated for the four mod- 
els from Table [T] are presented in Fig. [TJa). We see that 
all cumulative XLF are harder than the observed XLF at 
luminosities below 10 38 erg s _1 . This fact is easy to under- 
stand from our Fig. [5] and equation (|2ip : the typical slope 
of 1.25 is related to the intermediate characteristic ages (see 
cq.[7]), where the efficiency varies strongly. The position of 
the cutoff depends not only on the mean birth period, but 
also on the width of the distribution. For example, model 
3 has the largest mean period, but because of a large dis- 
persion, the XLF extends to very high luminosities without 
a visible break. On the other hand, model 2 has a rather 
small mean period, but the XLF cuts off sharply, because 
of the zero o p and the absence of fast pulsars. Model 1 has 
the largest number of bright pulsars because of the smallest 
mean period (see also Fig. [2}. Model 4 also shows a cutoff at 
rather small luminosity in spite of the small initial periods, 
because of the narrow magnetic field distribution. 

The number of high-luminosity pulsars depends on the 
maximum efficiency 770- Decreasing 770 leads to a smaller cut- 
off luminosity (compare thick and thin curves in Fig. [Zt a ))- 
However, if the cutoff is at very large luminosity (as e.g. in 
the case of models 1 and 3), variations in 770 do not affect 
significantly the observed XLF, at least in the range of lumi- 
nosities log L < 40. The XLF normalization scales linearly 
with the beaming factor b. 

Dividing the pulsar cumulative XLF by the observed 
XLF, we obtain the fraction of pulsars as a function of lumi- 
nosity. It is an increasing function of luminosity and reaches 
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Figure 8. Allowed region for the parameters of the birth period distribution derived from the XLF of sources in the nearby galaxies and 
the fraction of pulsars in the observed cumulative XLF jMineo et alj|2012h . Parameters along the bold red lines satisfy equation 1M6II . 
The solid, dashed, dotted and dot-dashed lines correspond to different pairs of (b,r]o) = (1,1), (1,0.3), (0.3,1) and (0.3,0.3), respectively. 
Regions to the right of these lines are allowed and to the left are forbidden. Dashed black and dotted blue contours give the pulsar 
fraction for b = 1 at luminosities above 10 39 erg s — 1 and 10 40 erg s — 1 , respectively. Calculations are performed for the average values 
log(iJ) = 12.6, Ob = 0.4 and K = 0.01. Positions of parameters listed in Table[T]are marked by different symbols. 



the maximum at log L between 38 and 41, depending on the 
model parameters. The maximum pulsar fraction reaches 
(0.2-0.5)6 for all models. 



the range of periods from tenths t o hundreds of milliseconds. 
Only the distribution found by iFaucher-Giguere fc KaspH 
(2006) satisfies this condition. 



3.3.3 Constraints on the birth period distribution 

The differential XLF of the pulsars should not exceed the 
observed XLF at any luminosity. This condition gives us the 
opportunity to find constraints on the birth period distribu- 
tion. We can find the dependence a p ({po)), which satisfies 
the condition 



max 

L<L m „ 



XLF pu i sars (L, (logpo), o-p) 



XLF„ 



= 1, 



(46) 



where L max = 10 41 erg s 1 corresponds t o the maximal lu- 
minosity considered bv lMineo et al.l(|2012r i. As with the con- 
straints from the SNe, here the results depend on the as- 
sumed values of rjo and b (see Fig. [HJ. As we see, the mean 
birth period (for 6=1 and 770 = 1) has to be larger than 
25-250 ms, depending on the w idth of the period d is tribu - 
tion . Parameters consid ered by lArzoumanian et all (|2002l ) 
and lTakata et all (|201ll ) lie in the forbidden region. Parame- 
ters from other papers listed in Table [J are in the allowed re - 
gion. However, as it was shown bv lPopov fc Turollal (|2012h . 
the period distribution has to be rather wide and to cover 



For a smaller value of the maximal efficiency rjo = 0.3, 
the critical line (red bold dashed line in Fig. [8} shifts to the 
left and depends weaker on a v . Thus, the decrease of the 
efficiency leads to shorter allowed periods. Variations in the 
beaming factor lead to a stronger effect. For b = 0.3 and 770 = 
1, the allowe d regio n extends beyond the parameters from 
iTakata et al.l ll20Tll), but still cannot reach the parameters 



from 



Arzoumanian et al.l (|2002T ). 



We note here that for simulations we used the aver- 
age value of the pulsar birth rate re = .01, while it is more 
than 10 and 2 times smaller in models of|X rzoumanian et al.l 
( 2002) and ITakata et alJ (|201lh . respectively. Thus all con- 
sidered models are in principle allowed if one corrects for 
different re. However, for small b and rjo the constraints com- 
ing from the SNe (see Fig. |Bp) are actually stronger and rule 
out model 1, with model 4 being only marginally consistent 
with the data. 
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Figure 9. The fraction of the pulsars in the observed cumulative XLF llMineo et al.ll2012h at luminosities above 10 39 erg s — 1 (panel a) 
and above 10 4t) erg s — 1 (panel b) (see Section l4,l|l as a function of a p . Solid, dotted, dashed and dot-dashed lines correspond to the mean 
period of log(po) = —1.8, —1.5, —1.0 and 0.5, respectively. The upper black and the lower blue lines give the pulsar fraction for 770 = 1 
and 0.3, respectively. Calculations are performed for the values log(S) = 12.6, ctb = 0.4 and k = 0.01. 



4 PULSAR CONTRIBUTION TO ULX 

4.1 Dependence on the birth period distribution 

The total number of the luminous pulsars with luminosities 
greater that 10 39 erg s" 1 is very similar in all four models 
1-4 (see Fig. Eta)): 

iV obs (log L > 39) w 0.3 x 6 x SFR [m q yr -1 ] . (47) 

At larger luminosities this number depends strongly on the 
period distribution and the maximum efficiency. For exam- 
ple, model 2 predicts less than 0.016 x SFR pulsars above 
10 40 erg s _1 because of the cutoff at ~ 2 x 10 39 erg s _1 in the 
XLF (see Fig.Et a ))- On the other hand, model 3 gives about 
0.086 x SFR very bright pulsars. The pulsar fraction in the 
observed XLF at logL > 39 can be as high as (0.2-0.3)6 for 
all models. At even higher luminosities, this fraction drops 
in models 2 and 4 and increases in models 1 and 4. 

We can also calculate the pulsar fraction dependence 
on the combination (log(po), cr p ). This fraction calculated for 
sources with luminosities log L > 39 and > 40 is presented as 
contours in Fig. [5] We see that the models 2 and 3 lie nearly 
on the same curve. The explicit dependence of the pulsar 
fraction on the parameters of the birth period distribution 
is shown in Fig. EI We see that for small initial mean periods, 
the pulsar fraction is nearly independent of a p because the 
XLF cuts off at very high luminosities. The situation changes 
dramatically at large (po): the narrow period distribution 
now predicts cutoff at low luminosity and the pulsar fraction 
is negligible. At large widths a p ~ 1, the pulsar fraction is 
still large because of the large extent of the XLF. Situation 
is similar for the cutoff luminosity logL = 40, but now for 
small <T P the XLF cuts off close to the limiting luminosity 
even for rather short initial periods and the pulsar fraction 
is small in that case. For large a p ~ 1, the pulsar fraction 
exceeds that fraction for log L — 39 if log(po) < —0.5. 



4.2 Dependence on the maximum efficiency 

In a general case, the dependence of the pulsar fraction 
above logL = 39 on 770 can be easily seen in Fig. Ht a ). For 
small initial periods, the pulsar fraction is nearly indepen- 
dent of 7/0 for all a p because the XLF extends to very high 
luminosities. For larger initial mean periods, the dependence 
on 770 is strong for narrow distributions a p ~ 0, because the 
XLF has a sharp cutoff around log L = 39. For broad initial 
distributions with a p ~ 1, the pulsar fraction is still rather 
large and the dependence on 770 is not so strong. 

The pulsar fraction at luminosities in excess of log L = 
40 (see Fig.^b)) shows a similar behaviour, but dependence 
on 770 is stronger because typically the XLF cuts off at that 
luminosity even for small initial periods and large a p . 

The dependence of the pulsar fraction on the maximum 
efficiency 770 for the models listed in Table [T] is presented in 
Fig. 1101 Most of the models predict rather flat dependence 
on 770 of the pulsar fraction at logL > 39 because of the 
wide initial period distribution producing the XLF extend- 
in g to rather high lumi nosities. The only exception is model 
2 (|Gonthier et al.|[2002j 'l , which predicts a significant drop in 
the pulsar function below 770 = 0.5. This is a direct conse- 
quence of the fact that this model has a narrow period distri- 
bution and its XLF has a sharp cutoff at about logL = 39.3 
for 770 = 1.0. Thus we see that for a rather wide range of 770 
between 0.3 and 1 the pulsar fraction above logL = 39 is 
between about 10 and 30 per cent for all models. Obviously, 
the beaming can reduce this fraction proportionally and for 
6 — 0.3 it is then at least 3 per cent. 

For models 1 and 3, the pulsar fraction is even larger at 
very high luminosities in excess of logL = 40 reaching 0.46 
and 0.36, respectively. In those cases, the dependence on 770 
is also not strong. While for model 2, the pulsar fraction is 
below 3 per cent and scales approximately as 770- From Fig. [5] 
it is clear that the closer parameters of the birth periods are 
to the limiting (bold red) line, the larger is the pulsar frac- 
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Figure 10. The predicted fraction of pulsars at luminosities 
above 10 39 erg s — 1 (thick lines) and above 10 40 erg s — 1 (thin 
lines) in the aver age XLF of sources in the nearby galaxies 
llMineo et alj|2012ft taken in the form H43H as a function of the 
maximum efficiency. Black solid, blue dotted, red dot-dashed and 
green dashed curves correspond to models 1-4 from Table [T] re- 
spectively. 



tion. For large <r p , the cumulative XLF is less steep than the 
observed XLF and therefore the pulsar fraction is a mono- 
tonically growing function that can reach 100 per cent above 
log L = 40. Thus it is possible that the pulsar fraction among 
the brightest ULX is significantly larger than 10 per cent. 
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Figure 11. Distribution of pulsars over magnetic field and ini- 
tial periods. Black, red, green and blue solid curves encircle 90 
per cent of pulsars for models 1—4 from Table [T] respectively. The 
dotted and dashed curves encircle 90 per cent of initial pulsar dis- 
tribution that can be observed to radiate above 10 39 erg s — 1 and 
10 40 erg s — 1 , respectively. The mean and the standard deviation 
describing these distributions are given in Table [4] Dotted pink 
curves give the dependence r c = ri and r c = T2- Solid brown 
curves are the lines of constant X-ray luminosity of 10 39 and 
10 40 erg s — 1 for 770 = !■ 



4.3 Distribution functions of the luminous pulsars 

In order to describe the possible observational appearance 
of the pulsars that can be observed as ULXs, we find a pos- 
teriori distribution of pulsars with luminosities log L > 39 
and 40 over magnetic field and birth periods as well as over 
their true ages. Because the pulsar luminosity drops with 
time, we are interested only in pulsars that emit above a 
given limiting L at birth. On the logB-logpo plane, these 
are the pulsars to the left of the corresponding L — const 
line (see Figs [T] and [TTJ. The probability that a pulsar will be 
observed above a given luminosity threshold is proportional 
to the pulsar age when it crosses the limiting L = const 
line. Thus the density distribution of such pulsars on the 
logB-logpo plane (limited to the region left of the limit- 
ing luminosity line) is given by the product of the density 
distribution at birth and the true age t = (pi — po)/aB 2 : 

P(logpo,logB) oc (x H(p c -p ) x G(logpo;log(po),o-p) 
x G(logB;log(B),a s ), (48) 

where H is the Heaviside step function and p c is the pulsar 
period when it crosses the limiting luminosity line given by 
equation (flQ)) . 

These density distributions of the bright observed pul- 
sars for various models from Table Q] are shown in Fig. 1111 
These distributions are generally narrower than the orig- 
inal distribution and skewed towards smaller periods and 



larger magnetic fields. They are elongated along the line of 
constant luminosity. The mean values and the standard de- 
viations of these distributions are given in Table [4] We also 
find the distribution of true ages of these pulsars, which is a 
monotonically decreasing function and can be described by 
the median age r. 

The evolution of the luminosity for the average lumi- 
nous pulsar can be described as follows: during the first 
ten years the luminosity is nearly constant at the level of 
~ 10 40 erg s _1 . After that it starts to decrease and still ex- 
ceeds ~ 10 39 erg s _1 for the next 100 years, during which 
the pulsars can be observed as ULXs. After that the lumi- 
nosity decreases down to ~ 10 36 erg s _1 in about 1000 yr. 

The rotation-powered pulsars are often assumed to be 
non-variable sources, but there may exist some variability 
on the time scales shorter that ~ 100 yr related to the in- 
teraction of the SN remnant with t he PWN and the sur- 
round ing media. As it was shown bv lDwarkadas fc Gruszkol 
(2012), the SN remnant could show a variability at least on 
the time scales ~ 10 yr. This variability could depend on the 
scale and the spatial spectrum of inhomogeneities of the sur- 
rounding media, and the characteristic variability timescale 
may increase with the pulsar age. The spectrum in the 0.1- 
10 keV range should consist of the soft thermal component 
related to the shock and the power-law tail related to the 
synchrotron radiation both from the shocks and the central 
pulsar. 
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Table 4. Parameters of the birth period and the magnetic field 
distributions as well as the median true age of pulsars with ob- 
served luminosities in excess of a given value for models from 
Table [1] 



Model 


log(P0> 




log(B> 


0"B 


T 

yr 






log L > 39 






1 


-2.43 


0.24 


12.35 


0.34 


343 


2 


-1.52 


0.0 


13.13 


0.16 


96 


3 


-2.05 


0.39 


12.79 


0.44 


136 


1 


-1.80 


0.08 


12.67 


0.09 


105 






logL 


> 40 






1 


-2.46 


0.23 


12.39 


0.32 


174 


2 


-1.52 


0.0 


13.37 


0.13 


13 


3 


-2.22 


0.34 


12.71 


0.41 


52 


1 


-1.86 


0.07 


12.70 


0.08 


36 



5 SUMMARY 

In the present paper we have investigated the question 
whether rotation-powered pulsars and PWN could be ob- 
served as some subclass of ULXs, and, if it is so, what is the 
fraction of pulsars in the whole ULX population. 

First, we developed an analytical model of the X-ray 
luminosity function, by solving the evolution equation for 
the period distribution of the pulsars. We derived both the 
steady state and the time-dependent solution. The steady- 
state solution is transformed to the pulsar XLF. We showed 
that this XLF has a broken power-law shape, reflecting the 
complex behaviour of the efficiency, with the high luminos- 
ity cut-off which location and shape are determined by the 
parameters of the birth period and magnetic field distribu- 
tions. The location of the cutoff mainly depends on the mean 
birth period. For short enough birth periods, the cutoff may 
lie above 10 39 erg s _1 . Therefore, the existence of luminous 
pulsars is possible. 

The time-dependent solution tells us about the evolu- 
tion of the distribution functions of the pulsars. We have 
shown that at large ages the period distribution becomes 
a delta- function- like peaking at \fa±B. The time-dependent 
luminosity distribution is more complicated due to the com- 
plexity of the luminosity-period relation. It can be multi- 
modal with different modes related to the different regimes 
of the efficiency of conversion of the rotation energy losses 
to the X-ray radiation. As the age of the pulsars increases, 
the luminosity distribution becomes more symmetric. 

We found constraints on the parameters of the birth 
period distribution using the observed XLF of the source s 
in the nearby galaxies obtained by iMineo et al.l l|2012f ). 
We found that the mean birth period cannot be shorter 
than 10-30 ms, depending on the width of the distribution. 
There fore, the parameters derived by lArzoumanian et al.l 
(2002) lie in the forbidden region for the typically assumed 
pulsar production rates . Acco unting for the recent find- 
ings oflPopov fc Turollal (120121) the parameters obtained by 



iFaucher-Giguere fc Kaspl (|2006t ) are the most reliable. 

We discussed the influence of the beaming and the max- 
imal efficiency on the luminosity function. For our calcula- 
tions we assumed conservatively b > 0.3, but the results can 
be easily scaled to the different values of the beaming. The 



number of the observed pulsars and their contribution to the 
ULX population depend linearly on b. The influence of the 
maximal efficiency is more complex, because it affects only 
the high luminosity tail of the pulsar XLF. For 770 = 0.3 
the allowed parameter space of the birth period distribution 
expands towards the shorter bir th periods. The fra ction of 
pulsars in the observed XLF of IMineo et all (|2012h would 
be smaller for smaller values of the efficiency and it strongly 
depends on the luminosity above which this fraction is com- 
puted. We showed that for broad initial period distributions, 
the pulsar fraction is a weak function of 770. 

We have also obtained constraints on the period dis- 
tribut ion by applying the method proposed by iPerna et al.l 
|2008). We derived the luminosity function of core-collapse 
SNe, using published X-ray light curves and compared it 
to the time-dependent luminosity function for pulsars. We 
found that the observed luminosities of the SNe are consis- 
tent with the mean birth period of po >0. 015-1 s, depend- 
ing on the width of the distribution, maximum efficiency 
and the beaming fact or. These con s traint s are in agreement 
with those derived bv lPerna et al. I (HH). 

We estimated a possible fraction of the pulsars in the 
whole population of ULX, using the observed XLF from 
M mco et al.li|2012T ). For the models considered in the previ- 
ous studies of pulsar populations, the predicted fraction of 
luminous pulsars can be in excess of 3 per cent for the sources 
with luminosities greater than 10 39 erg s -1 . At this moment , 
about 500 ULXs have been discovered l| Walton et al.ll201ll : 
ISwartz et al.l [20T1I ; iFeng fc Soriall201ll ) and we expect that 
at least ~ 15 of those should be associated with the rotation- 
powered pulsars. The models predict the pulsar fraction 
above 10 40 erg s _1 at the level of 1-40 per cent. 

Therefore, we might potentially observe bright pulsars 
as ULXs in galaxies with high SFR. These pulsars should 
have almost constant luminosity during the first hundred 
years after their birth, but there may exist some variability 
on the timescale of ~10 yr related to the interaction of the 
expanding SN remnant shell with the surrounding media. 
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